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In this article, we proposed a susceptible-infected model with identical infectivity, in which, at 
every time step, each node can only contact a constant number of neighbors. We implemented this 
model on scale-free networks, and found that the infected population grows in an exponential form 
with the time scale proportional to the spreading rate. Further more, by numerical simulation, we 
demonstrated that the targeted immunization of the present model is much less efficient than that 
of the standard susceptible-infected model. Finally, we investigated a fast spreading strategy when 
only local information is available. Different from the extensively studied path finding strategy, 
the strategy preferring small-degree nodes is more efficient than that preferring large-degree nodes. 
Our results indicate the existence of an essential relationship between network traffic and network 
epidemic on scale-free networks. 
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I. INTRODUCTION 



Since the seminal works on the small-world phe- 
nomenon by Watts and Strogatz [J] and the scale-free 
property by Barabasi and Albert [3 , the studies of com- 
plex networks have attracted a lot of interests within the 
physics community y, ll|- One of the ultimate goals 
of the current studies on complex networks is to un- 
derstand and explain the workings of the systems built 
upon them ja, |6j . The previous works about epidemic 
spreading in scale-free networks present us with com- 
pletely new epidemic propagation scenarios that a highly 
heterogeneous structure will lead to the absence of any 
epidemic threshold (see the review papers |3, |^ and 
the references therein). These works mainly concen- 
trate on the susceptible-infected-susceptible (SI S) |9lll(]| 
and susccptible-infected-removed (SIR) jT^, [l2 models. 
However, many real epidemic processes can not be prop- 
erly described by the above two models. For example, in 
many technological communication networks, each node 
not only acts as a communication source and sink, but 
also forwards information to others |13l [lj| ■ In the pro- 
cess of broadcasting [lallaj each node can be in two dis- 
crete states, either received or unreceived. A node in the 
received state has received information and can forward 
it to others like the infected individual in the epidemic 
process, while a node in the unreceived state is similar 
to the susceptible one. Since the node in the received 
state generally will not lose information, the so-called 
susceptible-infected (SI) model is more suitable for de- 
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scribing the above dynamical process. Another typical 
situation where the SI model is more appropriate than 
SIS and SIR models is the investigation of the dynamical 
behaviors in the very early stage of epidemic outbreaks 
when the effects of recovery and death can be ignored. 
The behaviors of the SI model are not only of theoreti- 
cal interest, but also of practical significance beyond the 
physics community. However, this has not been carefully 
investigated thus far. 

Very recently, Barthelemy et al. [i3ill3 studied the SI 
model in Barabasi- Albert (BA) scale- free networks Q, 
and found that the density of infected nodes, denoted 
by i{t), grows approximately in the exponential form, 
i{t) ~ e'^*, where the time scale c is proportional to 
the ratio between the second and the first moments of 
the degree distribution, c ^ {k'^)/{k) . Since the degree 
distribution of the BA model obeys the power-law form 
P{k) ~ k~^ with 7 = 3, this epidemic process has an infi- 
nite spreading velocity in the limit of infinite population. 
Following a similar process on random Apollonian net- 
works il9i, I20I, i2l| and the Barrat-Barthelemy-Vespignani 
networks ;22t t23j ■ Zhou et al. investigated the effects 
of clustering [13 and weight distribution [23 on SI epi- 
demics. And by using the theory of branching processes, 
Vazquez obtained a more accurate solution of i{t), in- 
cluding the behaviors with large t "2^. The common as- 
sumption in all the aforementioned works [13, HM Il9l |24J 
is that each node's potential infection-activity (infectiv- 
ity), measured by its possibly maximal contribution to 
the propagation process within one time step, is strictly 
equal to its degree. Actually, only the contacts between 
susceptible and infected nodes have possible contribu- 
tions in epidemic processes. However, since in a real epi- 
demic process, an infected node usually does not know 
whether its neighbors are infected, the standard network 



SI model assumes that each infected node will contact 
every neighbor once within one time step p?7|. thus the 
infectivity is equal to the node degree. 

The node with very large degree is called a hub in 
network science y, |j, |M ly|> while the node with great 
infectivity in an epidemic contact network is named su- 
perspreader in the epidemiological literature J2a . l27l |28| . 
All the previous studies on SI network model have a basic 
assumption, that is, hub = super spreader. This assump- 
tion is valid in some cases where the hub node is much 
more powerful than the others. However, there are still 
many real spreading processes, which can not be properly 
described by this assumption. Some typical examples are 
as follows. 

• In the broadcasting process, the forwarding capacity 
of each node is limited. Especially, in wireless multihop 
ad hoc networks, each node usually has the same power 
thus almost the same forwarding capacity |29l |. 

• In epidemic contact networks, the hub node has many 
acquaintances; however, he/she could not contact all 
his/her acquaintances within one time step. Analogously, 
although a few individuals have hundreds of sexual part- 
ners, their sexual activities are not far b eyon d a normal 
level due to the physiological limitations [Sfl IsJ |3l| . 

• In some email service systems, such as the Gmail sys- 
tem schemed out by Google "3^ , one can be a client only 
if he/she received at least one invitation from some ex- 
isting clients. And after he/she becomes a client, he/she 
will have the ability to invite others. However, the max- 
imal number of invitations he/she can send per a certain 
period of time is limited. 

• In network marketing processes, the referral of a 
product to potential consumers costs money and time 
(e.g. a salesman has to make phone calls to persuade his 
social surrounding to buy the product). Thus, generally 
speaking, the salesman will not make referrals to all his 
acquaintances 34|. 

In addition, since the infectivity of each node is as- 
signed to be equal to its degree, one cannot be sure which 
(the power-law degree distribution, the power-law infec- 
tivity distribution, or both) is the main reason that leads 
to the virtually infinite propagation velocity of the infec- 
tion. 



II. MODEL 

Different from the previous works, here we investigate 
the SI process on scale- free networks with identical infec- 
tivity. In our model, individuals can be in two discrete 
states, either susceptible or infected. The total popula- 
tion (i.e. the network size) N is assumed to be constant; 
thus, if S{t) and I{t) arc the numbers of susceptible and 
infected individuals at time t, respectively, then 



N = S{t)+I{t). 
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FIG. 1: (Color online) The infected density i{t) vs time, where 
i{t) — I{t)/N. The black and red curves result from the stan- 
dard SI network model and the present model. The numerical 
simulations are implemented based on the BA network '2] of 
size A'^ = 5000 and with average degree (k) = 6. The spread- 
ing rate is given as A = 0.01, and the data are averaged over 
5000 independent runs. 



during one time step. Accordingly, one can easily obtain 
the probability that a susceptible individual x will be 
infected at time step t to be 



A,(t) = l-(1-A)''(-^*- 



1) 



(2) 



where 6{x,t—l) denotes the number of contacts between 
X and the infected individuals at time t — I. For small A, 
one has 



X^{t)^X9{x,t-l). 



(3) 



Denote by A the spreading rate at which each susceptible 
individual acquires infection from an infected neighbor 



In the standard SI network model |17i ITM Il9l |. each 
infected individual will contact all its neighbors once at 
each time step, thus the infectivity of each node is de- 
fined by its degree and 9{x,t) is equal to the number of 
its infected neighbors at time t. In the present model, 
we assume every individual has the same infectivity A, 
in which, at every time step, each infected individual will 
generate A contacts where A is a constant. Multiple con- 
tacts to one neighbor are allowed, and contacts between 
two infected ones, although having no effect on the epi- 
demic dynamics, are also counted just like the standard 
SI model. The dynamical process starts by selecting one 
node randomly, assuming it is infected. 



III. SPREADING VELOCITY 

In the standard SI network model, the average infectiv- 
ity equals the average degree (fc). Therefore, in order to 
compare the proposed model with the standard one, we 
set A = (k). As shown in Fig. 1, the dynamical behav- 
iors of the present model and the standard one are clearly 
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FIG. 2: (Color online) The infected density i(t) vs time 
in normal (a) and single-log (b) plots. The black solid, 
red dot, green dash and blue dash-dot curves correspond to 
A = 0.01,0.001,0.0005 and 0.0001, respectively. In single- 
log plot (b), the early behavior of i{i) can be well fitted by 
a straight line, indicating the exponential growth of infected 
population. The inset shows the rescaled curves i{Xt). The 
four curves for different A collapse to one curve in the new 
scale Xt. The numerical simulations are implemented based 
on a BA network of size A'^ = 5000 and with average degree 
(k) = 6, and the data are averaged over 5000 independent 
runs. 



different: The velocity of the present model is much less 
than that of the standard model. 

In the following discussions, we focus on the proposed 
model. Without loss of generality, we set ^ = 1. Denote 
by ifc(i) the density of infected fc-degree nodes. Based on 
the mean-field approximation, one has 



d^fe(i) 



dt 



= Afcfl 



*^w]E^ 



k'Pik')tk'{t) 

k' Y,k" k"P{k") 



(4) 



where P{k) denotes the probability that a randomly se- 
lected node has degree k. The factor p- accounts for the 
probability that one of the infected neighbors of a node. 



FIG. 3: (Color online) The infected density i{t) vs time for 
different 7. The black squares, red circles, blue up-triangles, 
green down-triangles, and pink diamonds (from up to down) 
denote the cases of 7 = 2.0, 2.5, 3.0, 3.5 and 4.0, respectively. 
The numerical simulations are implemented based on the 
scale-free configuration network model. The networks are of 
size N = 1000 and with average degree (k) = 6, the spread- 
ing rate is given as A = 0.01, and the data are averaged over 
10000 independent runs. 



with degree k' , will contact this node at the present time 
step. Note that the infected density is given by 



iit) = j2^kit)Pik), 



so Eq. (4) can be rewritten as 
dife(t) Afc 



dt 



(fc) 



[l-ikit)]i{t). 



(5) 



(6) 



Manipulating the operator J^k -P(^) ^^ both sides, and 
neglecting terms of order O(i^), one obtains the evolution 
behavior of i{t) as follows: 



^{t) 



(7) 



where c oc A is a constant independent of the power-law 
exponent 7. 

In Fig. 2, we report the simulation results of the 
present model for different spreading rates ranging from 
0.0001 to 0.01. The curves i{t) vs t can be well fitted 
by a straight line in single-log plot for small t with slope 
proportional to A (see also the inset of Fig. 2b, where the 
curves for different values of A collapse to one curve in 
the time scale Xt), which strongly supports the analytical 
results. Furthermore, based on the scale-free configura- 
tion model [3^ yg , we investigated the effect of network 
structure on epidemic behaviors. Different from the stan- 
dard SI network model [i3, UM i which is highly affected 
by the power-law exponent 7, as shown in Fig. 3, the 
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FIG. 4: (Color online) The infected density i{t) va time with 
different vaccinating ranges. Figure 4a and 4b show the re- 
sults of targeted immunization for the standard SI process 
in normal and single-log plots, respectively. Correspondingly, 
figure 4c and 4d display the results for the present model. In 
all the four panels, the black solid, red dash, blue dot and 
green dash-dot curves represent the cases of / = 0, 0.001, 
0.005 and 0.01, respectively. The numerical simulations are 
implemented based on a BA network of size N = 5000 and 
with average degree (fc) = 6, the spreading rate is given as 
A — 0.01, and the data are averaged over 5000 independent 
runs. For comparison, the infectivity of the present model is 
set as ^ = (fc) = 6. 



exponent 7 here has almost no effects on the epidemic 
behaviors of the present model. In other words, in the 
present model, the spreading rate A, rather than the het- 
erogeneity of degree distribution, governs the epidemic 
behaviors. 



IV. TARGETED IMMUNIZATION 

An interesting and practical problem is "whether the 
epidemic propagation can be effectively controlled by vac- 
cination aiming at part of the population [3, H |33 . The 
most simple case is to select some nodes completely ran- 
domly, and then vaccinate them. By applying the perco- 
lation theory, this case can be exactly solved [3^ l39l| . 
The corresponding result shows that it is not an ef- 
ficient immunization strategy for highly heterogeneous 
networks such as scale-free networks. Recently, some ef- 
ficient immunization strategies for scale-free networks are 
proposed. On the one hand, if the degree of each node can 
not be known clearly, an efficient strategy is to vaccinate 
the random neighbors of some randomly selected nodes 
since the node with larger degree has greater chance to 
be chosen by this double-random chain than the one with 
small degree [iy, |41| . On the other hand, if the degree 
of each node is known, the most efficient immunization 
strategy is the so-called targeted immunization |42L |43J , 
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FIG. 5: (Color online) The infected density i{t) vs time for 
different /3. In Figure 5(a), the black solid, blue dot, magenta 
dash-dot, red dash and green dash-dot-dot curves correspond 
to /3 = 0, -1, -2, 1 and 2, respectively. In Figure 5(b), the black 
solid, red dash, blue dot, green dash-dot, magenta dash-dot- 
dot and cyan short-dash curves, from up to down, correspond 
to /3 = 0, 0.1, 0.2, 0.3, 0.4 and 0.5, respectively. In Figure 
5(c), the black solid, red dash, blue dot, green dash-dot, ma- 
genta dash-dot-dot and cyan short-dash curves correspond to 
/3 = 0, -0.1, -0.2, -0.3, -0.4 and -0.5, respectively. The numeri- 
cal simulations are implemented based on the extensional BA 
network of size A'^ = 5000 and with average degree (k) = 6, 
the spreading rate is given as A = 0.01 and the data are av- 
eraged over 5000 independent runs. 



wherein the nodes of highest degree are selected to be 
vaccinated (see also a similar method in Ref. 44]). 

Here, we compare the performance of the targeted im- 
munization for standard SI model and the present model. 
To implement this immunization strategy, a fraction of 
population having highest degree, denoted by /, are se- 
lected to be vaccinated. That is to say, these Nf nodes 
will never be infected but the contacts between them and 
the infected nodes are also counted. Clearly, in both the 
two models, the hub nodes have more chances to receive 
contacts from their infected neighbors, thus this targeted 
immunization strategy must slow down the spreading ve- 
locity. In Fig. 4a and Fig. 4b, we report the simulation 
results for the standard SI model. The spreading ve- 
locity remarkably decreases even only a small fraction, 
/ = 0.001, of population get vaccinated, which strongly 
indicate the efficiency of the targeted immunization. Rel- 
atively, the effect of the targeted immunization for the 
present model is much weaker (see Fig. 4c and Fig. 4d). 
The difference is more obvious in the single-log plot (see 
Fig. 4b and Fig. 4d): The slope of the curve i{t) ~ i, 
which denotes the time scale of the exponential term 
that governs the epidemic behaviors, sharply decreases 
even only a small amount of hub nodes are vaccinated in 
standard SI process while changes slightly in the present 
model. 



V. FAST SPREADING STRATEGY 

As mentioned in the Sec. 4, previous studies about 
network epidemic processes focus on how to control the 
epidemic spreading, especially for scale-free networks. 
Contrarily, few studies aim at accelerating the epidemic 
spreading process. However, a fast spreading strategy 
may be very useful for enhancing the efficiency of net- 
work broadcasting or for making profits from network 
marketing. In this section, we give a primary discussion 
on this issue by introducing and investigating a simple 
fast spreading strategy. Since the whole knowledge of 
network structure may be unavailable for large-scale net- 
works, here we assume only local information is available. 

In our strategy, at every time step, each infected node 
X will contact its neighbor y (in the broadcasting process, 
it means to forward a message to node y) at a probabil- 
ity proportional to fc^, where ky denotes the degree of y. 
There are two ingredients simultaneously affect the per- 
formance of the present strategy. On the one hand, the 
strategy preferring large-degree node (i.e. the strategy 
with /3 > 0) corresponds to shorter average distance in 
the path searching algorithm \43L |4^ , thus it may lead 
to faster spreading. On the other hand, to contact an al- 
ready infected node (i.e. to forward a message to a node 
having already received this message) has no effects on 
the spreading process, and the nodes with larger degrees 
are more easily to be infected according to Eq. (6) in the 
case of /3 = 0. Therefore, the strategy with /3 > will 
bring many redundant contacts that may slow down the 



spreading. For simplicity, we call the former the shorter 
path effect (SPE), and the latter the redundant contact 
effect (RCE). 

Figure 5(a) shows the density of infected individuals 
i{t) as a function of t for different /?. Clearly, due to the 
competition between the two ingredients, SPE and RCE, 
the strategies with too large (e.g. /3 = 1, 2) or too small 
(e.g. P = —1, —2) /3 are inefficient comparing with the 
unbiased one with f3 — 0. The cases when f3 is around 
zero are shown in Figs. 5(b) and 5(c). In Fig. 5(b), one 
can see that the RCE plays the major role in determining 
the epidemic velocity when /3 > 0; that is, larger /3 leads 
to slower spreading. As shown in Fig. 5(c), the condition 
is much more complex when /3 < 0: In the early stage, 
the unbiased strategy seems better; however, as time goes 
on, it is exceeded by the others. 



VI. CONCLUSION AND DISCUSSION 

Almost all the previous studies about the SI model 
in scale-free networks essentially assume that the nodes 
of large degrees are not only dominant in topology, but 
also the superspreaders. However, not all the SI network 
processes can be appropriately described under this as- 
sumption. Typical examples include the network broad- 
casting process with a limited forwarding capacity, the 
epidemics of sexually transmitted diseases where all in- 
dividuals' sexual activities are pretty much the same due 
to the physiological limitations, the email service sys- 
tems with limited ability to accept new clients, the net- 
work marketing systems where the referral of products 
to potential consumers costs money and time, and so on. 
Inspired by these practical requirements, in this article 
we have studied the behaviors of susceptible-infected epi- 
demics on scale-free networks with identical infectivity. 
The infected population grows in an exponential form in 
the early stage. However, different from the standard 
SI network model, the epidemic behavior is not sensitive 
to the power-law exponent 7, but is governed only by 
the spreading rate A. Both the simulation and analytical 
results indicate that it is the heterogeneity of infectivi- 
ties, rather than the heterogeneity of degrees, governs the 
epidemic behaviors. Further more, we compare the per- 
formances of targeted immunization on the standard SI 
process and the present model. In this standard SI pro- 
cess, the spreading velocity decreases remarkably even 
only a slight fraction of population are vaccinated. How- 
ever, since the infectivity of the hub nodes in the present 
model is just equal to that od the small-degree node, the 
targeted immunization for the present model is much less 
efficient. 

We have also investigated a fast spreading strategy 
when only local information is available. Different from 
previous reports about some relative processes taking 
place on scale- free networks [43, |43 , we found that the 
strategy preferring small-degree nodes is more efficient 
than those preferring large nodes. This result indicates 



that the redundant contact effect is more important than 
the shorter path effect. This finding may be useful in 
practice. Very recently, some authors suggested using 
a quantity named saturation time to estimate the epi- 
demic efficiency [47.. 48j , which means the time when the 
infected density, i{t), firstly exceeds 0.9. Under this cri- 
terion, the optimal value of j3 leading to the shortest 
saturation time is -0.3. 

Some recent studies on network traffic dynamics show 
that the networks will have larger throughput if u sing 
routing strategies preferring small-degree nodes [43, llfl 
l51| . It is because this strategy can avoid possible conges- 
tion occurring at large-degree nodes. Although the quan- 
titative results are far different, there may exist some 
common features between network traffic and network 



epidemic. We believe that our work can further enlighten 
the readers on this interesting subject. 
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